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ABSTRACT 


The  Dynamic  Model  of  Modern  Military  Conflict  developed 
by  Dr.  Paul  Moose,  Naval  Postgraduate  School,  is  described 
by  its  system  of  differential  equations  followed  by  an  inves¬ 
tigation  of  its  behavior.  This  investigation  is  predicated 
by  an  analysis  of  the  model's  stability  about  equilibrium 
using  a  method  attributed  to  the  study  of  ecosystems.  The 
basis  for  this  analysis  is  the  formulation  and  subsequent 
evaluation  of  a  community  matrix  termed  the  conflict  matrix 
for  this  investigation. 

Specific  investigation  of  the  model's  behavior  is  observed 
varying  the  Command,  Control,  Communications  and  Intelligence 
(C^I)  enhancements  of  one  combat  force  against  that  of  an 
opposing  force  while  one  force  utilizes  misinformation  and 
deception  as  a  counter  C^  tactic  and  the  other  uses  physical 

j 

destruction  as  a  counter  C^  tactic. 

Potential  military  use  of  this  model  as  an  analytical 
tool  for  playing  "what  if?"  type  wargames  is  envisioned 
following  further  research  and  study  of  the  effects  on  model 
behavior  by  varying  other  parameters  not  specifically  addressed 
in  this  initial  investigation. 
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I .  INTRODUCTION 


A  military  conflict  has  a  variety  of  elements  that 
potentially  influence  its  eventual  outcome.  Commanders  of 
the  various  echelons  within  the  combat  environment  determined 
by  a  particular  military  conflict  are  faced  with  decisions 
that  require  the  fusion  of  information  from  countless  sources. 
These  sources  may  be  highly  sophisticated  sensors,  relaying 
timely  digital  information  concerning  enemy  force  concentra¬ 
tions  or  movements,  aircraft  or  missile  tracks  both  friendly 
and  nonfriendly,  or  perhaps  intercepted  communications  signals 
These  sources  could  also  be  the  professional  schools,  exercise 
participation  or  intelligence  documents  that  have  molded  the 
commander's  decision  process  over  his  military  career  whether 
it  be  extensive  or  minimal.  Regardless  of  the  case,  decisions 
are  made  with  some  reliance  on  or  rejection  of  prior  knowledge 
The  assimulation  of  this  prior  knowledge,  whether  relative  to 
the  situation  or  not  can  be  a  monumental  task  in  itself. 

The  sophistication  of  computer  driven  information  systems 
and  a  presumed  need  for  rapid  decision  making  either  by  the 
man  as  a  commander  or  the  machine  as  his  extension  has  caused 
a  proliferation  of  model  development  in  recent  years  [Ref.  1] . 
These  models  attempt  to  reflect  the  variety  of  elements  that 
potentially  influence  the  outcomes  of  military  engagements, 
manipulate  them  according  to  some  predetermined  scheme  or 


algorithm  2Uid  finally  output  some  prescription  or  calculated 
result.  The  output  of  such  models,  together  with  the  parti¬ 
cular  input  elements  and  scheme  of  manipulation  offer 
limitless  applications  either  for  direct  command  decision 
making,  training  or  running  countless  "what  if"  exercises 
relative  to  known  or  suspected  enemy  and  friendly  combat 
information  or  characteristics. 

A.  PURPOSE 

The  primary  purpose  of  this  thesis  is  to  initiate  inves¬ 
tigation  of  one  such  model  that  exhibits  the  potential  for 
use  in  a  variety  of  areas,  but  in  particular  in  looking  at 
the  "what  if"  questions  associated  with  command,  control  and 
communications  (C^)  and  covmter  C^.  The  specific  model  under 
investigation  is  the  Dynamic  Model  of  Modern  Military  Conflict 
proposed  by  Dr.  Paul  Moose,  Academic  Associate  for  the  Joint 
Command,  Control  and  Communications  curricula  at  the  Naval 
Postgraduate  School  [Ref.  2]. 

B .  SCOPE 

The  basic  description  of  the  model  is  given  in  Chapter  2. 
This  description  is  in  the  form  of  a  system  of  differential 
equations  similar  to  that  which  describes  Lanchester-type 
attrition  models,  however  it  includes  terms  relative  to  C^ 
and  counter  C^  for  opposing  forces.  In  this  respect,  informa¬ 
tion  possessed  by  opposing  forces  is  elevated  to  an  equivalent 
level  with  the  forces  themselves.  Furthermore,  the  description 
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includes  a  discussion  on  how  the  equations  of  this  model  are 
also  similar  to  evolution  equations  which  describe  multi¬ 
species  ecosystems.  Such  systems  include  natural  death  (loss) 
rates  and  replenishment  of  species  within  a  population  thus 
leading  to  a  mathematical  discussion  of  this  model's  behavior 
relative  to  an  established  equilibrivun. 

Following  this  description.  Chapter  3  looks  at  the 
stability  of  the  model  in  terms  of  a  "community  matrix", 
another  reference  to  ecosystems  which  incorporates  those 
parameters  which  tend  to  increase  or  reduce  species  within 
a  population  or  commvinity.  This  "commvinity  matrix"  is 
renamed  the  "conflict  matrix"  for  this  model  and  explored 
by  assigning  a  certain  set  of  asymmetrical  values  to  the 
elements  of  the  matrix  relative  to  hypothetical  opposing 
forces  in  a  military  conflict.  Specifically  the  investiga¬ 
tion  undertaken  in  this  chapter,  demonstrates  conditions 
where  opposing  forces  stabilize  at  a  set  equilibrivun  point 
or  diverge  from  this  point  in  unstable  fashion. 

Chapter  4  looks  at  more  specific  behaviors  relative  to 
the  stability  and  instability  demonstrated  in  Chapter  3. 

The  identical  initial  conditions  used  in  Chapter  3  are 
again  used  in  this  chapter  and  represent  arbitrary  units 
rather  than  specific  calculated  values.  These  values  were 
chosen  purely  to  demonstrate  the  general  behavior  of  the 
model  rather  them  replicate  any  specific  military  conflict. 
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Finally  effects  on  equilibrium  relative  to  the  initial 
conditions  used  during  stability  and  behavior  investigations 
are  explored  in  Chapter  5.  Specifically,  the  initial  condi¬ 
tions  are  aimed  at  isolating  the  effects  of  varying  the 

3 

command,  control,  communications  emd  intelligence  (C  I) 
enhancements  of  one  force  against  an  opposing  force's  C^I 
enhancements.  Utilizing  these  results,  further  demonstration 
of  the  model  behavior  is  then  achieved  by  varying  resource 
and  information  variables  associated  with  opposing  military 
forces. 

In  essence,  this  thesis  is  concerned  with  am  initial 
demonstration  that  the  Dynamic  Model  of  Modern  Military 
Conflict  can  perform  as  might  be  expected  rather  than  prove 
that  it  works  in  all  cases.  This  approach  opens  several 
avenues  to  continued  research  leading  to  a  more  complete 
demonstration  of  its  usefulness  in  evaluating  the  influence 
of  C^I  amd  counter  C^  in  modern  military  conflict. 
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II.  A  DYNAMIC  MODEL  OF  MODERN  MILITARY  CONFLICT 

A.  DEFINITION 

The  Dynamic  Model  of  Modern  Military  Conflict  incorporates 

the  effects  of  C^I  with  basic  attrition  coefficients  such  as 

those  found  in  Lanchester-type  models.  Specifically  the 

effects  of  the  above  are  described  in  terms  of  information 

(I  and  I  )  and  general  purpose  resources  (M  and  M  )  for 

X  y  X  y 

opposing  X  and  Y  combat  forces  [Ref.  3].  The  system  of 
differential  equations  which  characterize  this  model  are 
given  below  followed  by  sections  in  this  chapter  which 
briefly  touch  on  the  background  of  the  model  and  provide 
specific  definitions  of  the  terms  involved. 

(Eqn  2.1) 

d<I^  /dt 

<Eqn  2.2) 

d(M^/dt  =  -  S^«^]  - 

(Eqn  2.3) 

d(I  )/dt  =  I[-A-aI-YM]+CIM+  QI 

y  '  y*"  y  yx  X  'yx  x-*  yy  ye  y  y  ye 

(Eqn  2.4) 

d(M^/dt  .  My(-By  -  -  dyyVy  -  -  V“ye^  *  We 
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B .  BACKGROUND 


This  system  of  four  first  order  non-linear  differential 
equations  is  of  the  form  =  -F(S)  +  2  ^ 

equations  in  this  form  are  knovm  as  evolution  equations. 

S  and  Q  are  4x1  column  vectors  and  there  are  four  func- 
tional  relationships  which  are  considered  to  be  at  most 
quadratic.  Expanding  the  loss  function  term  gives  the 
following : 


a.  .S  . , 
J-D  3 


(Eqn  2.5) 


1,2, 3, 4 


i  .e. 


-  (S)=  -  quadrative  term  -  linear  term  (Eqn  2.6) 

This  system  as  it  turns  out  is  a  very  generalized  form  of  the 
Lotka-Volterra  equations  for  multi-species  ecosystems  [Ref. 

4:  37-38].  More  specifically  it  takes  on  the  form: 

(Eqn  2.7) 


i  =  1,2, 3, 4 


C.  EXPLANATION  OF  TERMS 

In  general  the  greek  letter  coefficients  are  quadratic, 
and  the  others  are  linear.  Following  is  the  definition  of 
the  terms  that  characterize  the  model; 

Ij^  -  Information  of  force  X 
I^g  -  (Scune  as  above  at  equilibrium) 

Mjj  -  General  purpose  resources  for  force  X 
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M. 


xe 


B. 


yx 


yx 


XX 


XX 


a. 


yx 


'yx 


yx 


-  (Saioe  as  above  at  equilibrium] 


-  Natural  death  (loss)  rate  for  Cl  resources 
of  force  X 


-  Natural  death  (loss)  rate  for  general  purpose 
resources  of  force  X 


-  Lanchester  "area  fire"  coefficient,  X  on  Y 


-  Lanchester  "aimed  fire"  coefficient,  X  on  Y 


-  Fixed  replenishment  rate  for  C  I  resources  of 
force  X 


-  Fixed  replenishment  rate  for  general  purpose 
resources  of  force  X 


-  Replenishment  from  general  purpose  resources 
for  C3i  production,  force  X 


-  Diminishment  of  X  force  for  Cl  production 


-  Counter  C'^  by  force  X  on  force  Y  by  misinformation 
and  deception  from  c3i  resources 


-  Counter  C  by  force  X  on  force  Y  by  physical 
destruction  from  general  purpose  resources 


-  Cl  (quadratic)  force  effectiveness  enhancements, 
X  on  Y 


yx 


-  C^I  (linear)  force  effectiveness  enhancements, 
X  on  Y 


The  remaining  coefficients  with  y  subscripts  or  xy  subscripts 
refer  to  the  Y  force  or  Y  on  X  activities. 


D.  EXPECTED  MODEL  BEHAVIOR 


The  actual  behavior  of  the  model  is  dependent  initially 
on  the  equilibrium  for  the  system.  One  equilibrium  point  is 
determined  by  solving  the  following  linear  set  of  equations 


(Eqn  2.8) 


“q 

^xx 

-a 

xy 

_ 1 

^  m 

I 

xe 

V 

V 

0 

-<^xy^xy> 

-(V^xy> 

M 

xe 

- 

®x 

+ 

R 

X 

“^yx 

0 

'^yy 

I 

ye 

A 

y 

yx^ 

-(^yx^yx^ 

“^yy 

Q 

ye 

B 

. 

Once  equilibrixim  is  established  the  effects  of  varying 
certain  coefficients  on  dependent  variables  1^, 

can  be  illustrated  by  solving  the  system  of  differential 
equations  over  the  independent  variable  time. 


III.  STABILITY  ANALYSIS 


A.  DEFINITION 

Stability,  as  it  relates  to  a  deterministic  model  or 
deterministic  system  of  equations,  is  more  specifically 
defined  as  "neighborhood  stability"  or  "stability  in  the 
vicinity  of  equilibirum" .  For  population  models  with 
environmental  pareimeters  all  well  defined  constants, 
interest  centers  on  the  community  equilibrium  where  all 
species'  populations  have  time-independent  values,  i.e. 
all  net  growth  (death)  rates  are  zero.  Such  2m  equilibrium 
may  be  called  stable  if,  when  the  populations  are  perturbed, 
they  in  time  return  to  their  equilibrium  values;  the  return 
may  be  achieved  either  as  damped  oscillations  or  monotonically 
Conversely,  if  such  a  disturbance  tends  to  amplify  itself,  the 
system  is  called  unstable;  such  instability  may  appear  as  an 
oscillatory  or  as  a  monotonic  growth  in  the  disturbance. 

This  definition  is  sufficient  for  linearized  systems,  however, 
it  may  be  misleading  for  non-linear  systems  when  full  global 
stability  of  a  system  is  of  concern  [Ref.  4:15], 

Full  global  stability  is  said  to  be  characterized  by 
neighborhood  analysis  if  the  existence  of  a  Lyapunov  function 
can  be  determined  [Ref.  4:15].  The  existence  of  such  a 
fvmction  is  not  explored  in  this  thesis,  however,  regardless 
of  whether  such  a  function  exists  or  not,  neighborhood  analysi 
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is  within  the  scope  of  this  thesis  and  founded  in  the 
generation  of  a  conttnunity  matrix. 

B .  COMMUNITY  MATRIX 

The  community  matrix  is  said  to  both  summarize  the  system 
(its  elements  being  determined  by  the  interactions  between  and 
within  species  near  equilibrium)  and  sets  the  neighborhood 
stability  by  the  sign  of  its  eigenvalues  (all  negative  real 
parts  for  all  eigenvalues  imply  stability) .  Specifically  it 
is  generated  by  the  following  mathematical  derivation  [Ref.  4; 
19-22]  . 

Since  the  model  is  of  the  form 


d(S) 

-jI - !‘!>  ♦  S 


its  multispecies  population  dynamics  are  given  by  the  set 
of  four  equations ; 

(Eqn  3.1) 


d(S.  (t)) 

- 1- -  =  F^(S^(t),  S2(t),  S3(t),  S4(t)),  i  =  1,2, 3, 4 

The  growth  rate  of  the  i'"  species  at  time  t  is  given  by 

the  nonlinear  function  of  all  relevant  interacting 

* 

populations.  Equilibrium  populations,  S^  ,  follow  from 
the  4  algebraic  equations  obtained  by  setting  all  growth 
rates  at  zero: 


0 


F,(Si  , 


(Eqn  3.2) 
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(Eqn  3.3) 


For  each  population,  let 
S^tt)  *  S^*  +  s^(,n 
where  measures  the  initially  small  perturbations  to 

the  i^  population.  Expanding  each  of  the  basic  equations 
around  this  equilibrium  in  a  Taylor  series  and  discarding  all 
second  order  and  higher  terms  in  the  population  perturbations 
s,  a  linearized  approximation  is  obtained: 

d  s.(t)  4 

—A -  =  I  a..  S.(t)  CEqn  3.4) 

j=l  ^ 

Equivalently,  in  matrix  notation 


d  X(t)  _ 


=  A  XCt). 


(Eqn  3.5) 


where  X  is  the  4x1  column  Matrix  of  the  s^,  and  A  is 
the  four  by  four  "community  matrix"  whose  elements  a^^ 
describe  the  effect  of  species  j  upon  species  i  near 
equilibrium.  The  elements  a. .  depend  both  on  the  details 
of  the  original  equations  and  on  the  values  of  the  equilibrium 
population  through 

3F.  * 


The  partial  derivatives,  denoting  the  derivatives  of  F^ 
keeping  all  populations  except  Sj  constant,  are  then 
evaluated  with  all  populations  at  their  equilibrium  values. 
The  corresponding  community  matrix  for  the  Dynamic  Model  of 
Modern  Military  Conflict  is  called  the  "conflict  matrix" 
and  its  elements  are  given  in  the  following  matrix: 
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V  *  V 


XX 


XX 


yx 


X  xy 


-Y, 


-C6  -Hi  ) 
yx  yx 


yx 

-(3  +b  ) 
yx  yx 


-^\^yx^V 


yy 


-V 

-Vxy> 


yy 

-(B  +6  +6  ) 

y  yx  yx 


C.  METHOD  OF  ANALYSIS 

To  investigate  the  stadaility  of  the  model,  eigenvalues 
are  computed  relative  to  the  conflict  matrix  given  above  and 
a  set  of  initial  values  for  the  coefficients  that  comprise 
the  elements  of  the  matrix.  Two  cases  are  investigated; 
one  in  which  element  d^^^  is  held  constant  and  element  d^^^ 
is  varied  across  a  range  of  values  and  another  where  d 

xy 

is  held  constant  and  d  is  varied  across  a  range  of  values, 

yx 

The  particular  method  of  computing  eigenvalues  is  the  IMSL 
routine  EIGRF  [Ref.  5]. 


D.  INITIAL  CONDITIONS 

The  initial  conditions  for  investigating  the  two  cases 
are  given  in  Table  I .  These  initial  conditions  were  chosen 
to  illustrate  some  of  the  potential  behavior  that  this  model 
may  exhibit.  The  variable  parameters  are  C^I  enhancements 
for  both  X  amd  Y  forces.  Furthermore  the  model  is 
asymmetrical;  X  is  using  physical  destruction  by  forces  as 
a  counter  technique  and  Y  is  using  misinformation  and 

deception  as  a  counter  technique. 
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TABLE  I 


INITIAL  CONDITIONS  FOR  STABILITY  INVESTIGATION 


XX 


xy 


xy 


xy 


xy 


a 


yx 


xy 


'  yx 


xy 


yx 

Equilibrium 


"  ‘^xx  “  ‘^yy 

V  = 

^x  = 

byX  =  Q-Q 

l.QO  (misinformation  and  deception) 
Q.Q 


0.0 

1.0  (physical  destruction) 
variable  C^I  enhancements  (Y  on  X) 

3 

variable  C  I  enhancements  (X  on  Y) 
is  taken  to  be  the  vector  1  ,  i.e. 


I 

xe 


M. 


xe 


1 
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Substitution  of  the  initial  parameter  conditions  into 


the  conflict  matrix  gives  the  following  matrix: 


-1.5 

0 

0 

-V 


Q 

-1.5 

-l.Q 

-l.Q 


-1.0 

"‘^yx 

-1.5 


0 

-1.0 

0 

-1.5 


The  particular  cases  investigated  then  are : 

d  held 
yx 

by  0.05. 
d^y  held 
by  0.05. 


1.  d  held  constant  at  1  and  d  varied  from  0  to  2 
yx  xy 


2.  d  held  constant  at  1  and  d  varied  from  0  to  2 
xy  yx 


E.  RESULTS  AND  DISCUSSION 

The  eigenvalues  generated  by  case  1  are  presented  in 

Table  II.  It  can  be  seen  that  for  values  of  d  greater 

xy 

than  0.80  a  positive  real  eigenvalue  appears.  This  is 
indicative  of  instability.  Similarly  all  real  parts  of 
eigenvalues  for  d^^^  less  than  or  equal  0 . 80  are  negative 
thus  indicating  stability. 

The  eigenvalues  generated  by  case  2  are  presented  in 
Table  III.  For  values  of  d^^  greater  than  0.55  instability 
is  indicated  while  stability  is  indicated  for  values  of  d 

yx 

less  than  or  equal  to  0.55. 
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EIGENVALUES  FOR  VARIABLE  C^I  ENHANCEMENTS  (Y  ON  X) 


hmmhhhhmmhhmmmmmmhhmhmhhhmhhmhmmmhhhhmhhh 

ooooooooooooooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooooooooaooo 

ooooooooooooooooooooooooooooooooooooooooo 


9  ooooooooooooooooooooooooooooooooooooooooo 

z  p*m<Na>vo-»»“ooir><NoomrMa*tncMoovn»-p»«ovO(Nooif>»-\arMeo^ovors»r*<noMr>o\orM 
cq  <»o9M^voin9(Ni^ooor^vo^<*>rMo9«aov0tn»r>jf-o\aor«ir>tf (Mr>oeor>ir>9(N<-oaor^ 
KS  (•>p<»«N<M<N<NfMr'j<N<N^»-»“»«^^*“OOOOoooaio\ON<5\<7\0'0»wooeoeoaoooooeof*»p» 


U  rornrnrnpnm«»lrr(rorn<nfnrr)(*»r«><»^fn(»><»>r>v»1pr)rOro<SrMrMfSfMfM{N(N<N(NtN(N{NfStN<N(N 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  i  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


MMHHMMMHMMMMMMMMHHMMHHMMMMHMHMMHHMMHMHMMM 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

f*)  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

z  i^m<NO'vo»^ooinrMflo»r>o>io'iftCMooir»»“r*^OvorM<Mmw»®f'i\o03’eorof^^ir»o^oo 

O  ®®r«|<>»<M<N(N<NrM<N»“*-^^^<-^00000000000000»-»-^^^»“r-CMr«JCM 


H  ooooooo®oooooooooooooooooooooooooooooo®oo 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


HHHHHHMHHHHMMMMHHHHHHHHHMMMHMHMHHHHHHHHHH 

o»®<N®04f®®®fMr»<Nf^<Mr»<sr«»r>»®^OMn*-p»a’Or«>®or*«r>i®r^^(M^®p»g3 

®tnvn®®r‘r*r«*®®®®oo«-»-<N<N®®aj»in®®r*®®<y>oo»-(N<N®»in®®r»® 

<N  ooooooooooooooooooooooooooooooooooooooooo 

I  I  I  I  I  I  •  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  t  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 
Z  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
H  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

M  . . . . . . . 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  i  I  t  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


MMMMMMHMHMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM 

o.»®<N®o»®®®<Nr«-<Nr-<Nr^(Nr^<N®3’a'«r>^r*'3-or'f^op*3’(N®r~-JT<N»-’®r«“® 

ifiif)invo®r^r^r»'®®o>®oo»*^<NfN®«^3’3'»ri®®r*®®®oo^«N<N®wm®®r'® 


^  OoOqOOOOOOOOOOoOOOOOOOOOOOOoOoOoOOOOOOOOO 
Z  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

w  ooooooooooooooooooooooooooooooooooooooooo 

o  inin»rtvnvninir)®ir>tnvnifnrtin«r»»n®»r>»riif)ir  mvnminiriiriiniDinriininininiftiriininmir 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  (•  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


>4  ovnoino»novfio®o®oirto»noino«noiroinoinoinoinomoino»r»omoif)o 

<^4-00000000000000000000 


24 
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The  results  of  the  above  two  cases  are  sunuaarized 

graphically  in  Figures  1,  2,  and  3.  Figure  1  depicts  the 

generation  of  eigenvalues  with,  positive  real  parts 

(instability)  beyond  0.5  and  0.8  for  d  and  d 

yx  yx 

respectively,  and  negative  real  parts  elsewhere  (stability) . 
Figure  2  represents  the  other  real  eigenvalue  which  is 
always  negative  thus  not  contributing  to  instability.  Figure 
3  represents  the  positive  imaginary  part  of  the  remaining  two 
eigenvalues  which  turn  out  to  be  complex  conjugate  pairs  with 
real  part  equal  to  -1.5.  Existence  of  these  imaginary  parts 
indicate  that  the  system  does  in  fact  oscillate  while  dampen¬ 
ing  out. 

The  eigenvalues  generated  by  these  two  cases  lead  to 

more  specific  investigation  of  stability  and  instability. 

The  eigenvalues  have  located  the  regions  of  stability  and 

instability  so  now  by  using  the  appropriate  values  of  d^^ 

and  d  the  behavior  of  the  model  can  be  observed  relative 
yx 

to  perturbations  of  the  dependent  variables  over  time  for 
stable  and  unstable  cases. 


Real 


0  12  3  4 


Figure  2.  Eigenvalues  with  all  negative  real 
parts  (imaginary  part  =  0) . 
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Imaginary 


IV.  SPECIFIC  STABILITY  ANALYSIS 


A.  METHOD  OP  ANALYSIS 

More  specific  illustrations  of  stability  in  the  Dynauaic 
Model  of  Modern  Military  Conflict  are  achieved  by  perturbing 
the  dependent  variables  and  observing  the  behavior  of  the 
model  relative  to  time.  The  initial  values  used  in  the 
stcJaility  analysis  of  Chapter  3  are  substituted  into  the 
differential  equations  which  characterize  the  model.  These 
differential  equations  turn  out  to  be  the  following  for  the 
parameter  values  of  Table  I . 


X 

.  ss  1  ^  -  T 

dt 

•  V  aXi  , 

d  M.. 

A 

dt 

=  1.5  -  d, 

d  I  ^ 

_ 1. 

dt 

=  1.5  -  I, 

d  M. 

=  1.5  -  d 

xy  y  X  '  y  '  xy 


y  *  “-S) 


I  -  M  (M  +  0.5)  +  d 
yx  X  y  '  X  yx 


(Eqn  4.1) 
(Eqn  4.2) 
(Eqn  4.3) 

(Eqn  4.4) 


The  Interactive  Ordinary  Differential  Equations  package 

(lODE)  was  en^loyed  using  these  differential  equations  and 

representative  values  for  d  and  d  (stable  and 

xy  yx 

unstable  values) [Ref.  6].  Specifically,  plots  were  generated 
to  graphically  represent  two  behaviors: 


1.  Dependent  variables  against  the  independent  variable 

i.e.,I  ,M,  I  ,M,  against  time. 

'  X  X  y  y  ’ 

2.  Phase  plots  of  opposing  forces  i.e.,  M  against  M  . 
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B.  STABILITY  WITH  d  AT  a. 4 

xy 

Table  IV  shows  the  initial  values  for  the  dependent 

variables,  time  duration  6uid  increment  for  plotting  for 

each  of  fourteen  trial  runs  of  this  model.  Both  d  and 

xy 

d  were  held  constant  at  Q.4  and  l.Q  respectively.  In 
yx 

addition,  the  final  column  in  Table  IV  represents  the  time 

at  which  all  dependent  variables  damped  out  to  or  very  nearly 

to  the  established  equilibrium  of  1. 

Graphically,  Figure  4  depicts  the  general  plot  where 

only  one  dependent  variable  was  perturbed  (trials  1  through 

4) .  Figures  5  through  8  further  illustrate  phase  plots  for 

these  S2une  trials  of  M  vs  M  illustrating  return  to 

X  y 

equilibrium  over  time. 

Two  other  trials  are  illustrated  in  Figures  9  through  12. 

Specifically  Figures  9  and  IQ  represent  perturbations  of  both 

I  and  I  (trial  6) .  The  interesting  phase  plot  in  Figure  10 
X  y 

is  indicative  of  the  oscillatory  nature  of  both  M  and  M 

X  y 

relative  to  these  perturbations.  Figure  11  depicts  a  triple 

perturbation  involving  1^^,  1^  and  M^,  with  the  resulting 

oscillatory  nature  of  M  and  M  shown  in  Figure  12  (Trial  13). 

X  y 

Although  the  above  illustrations  are  not  all  inclusive,  it 
should  be  noted  that  all  dependent  variables  did  in  fact  damp 
out  over  the  set  time  period  and  phase  plots  of  M„  cuid 

X  y 

revealed  return  to  equilibrium  in  all  trials. 
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TABLE  IV 

INITIAL  CONDITIONS  WITH 

d_  AT  Q.4  AND  d  AT  1.0  (STABLE) 
xy  yx 


Trial  # 

Time 

Increment 

5c 

Time  for 
90%  Damp 
out  (sec). 

1 

0 

to 

50 

.05 

1.25 

1.00 

1.00 

1.00 

7 

2 

0 

to 

50 

.05 

1.00 

1.25 

1.00 

1.00 

12 

3 

0 

to 

50 

.05 

1.00 

1.00 

1.25 

1.00 

9 

4 

Q 

to 

50 

.05 

l.QO 

1.00 

1.00 

1.25 

10 

5 

0 

to 

50 

.05 

1.25 

1.25 

1.00 

1.00 

17 

6 

0 

to 

50 

.05 

1.25 

1.00 

1.25 

1.00 

1 

7 

0 

to 

50 

.05 

1.25 

1.00 

1.00 

1.25 

1.5 

8 

0 

to 

50 

.05 

l.QO 

1.25 

l.QO 

1.25 

1.5 

9 

0 

to 

50 

.05 

l.QO 

1.00 

1.25 

1.25 

15.5 

10 

0 

to 

50 

.05 

1.00 

1.25 

1.25 

1.00 

4.0 

11 

0 

to 

50 

.05 

1.25 

1.25 

1.25 

1.00 

12.0 

12 

0 

to 

50 

.05 

1.25 

1.25 

1.00 

1.25 

11.0 

13 

0 

to 

50 

.05 

1,25 

1.00 

1.25 

1.25 

12.0 

05 


1.00  1.25  1.25  1.25  4.0 


Figure  4.  Single  Perturbation,  stable  (general  plot). 


.25) 


Figure  7.  Single  Perturbation,  stable  phase  plot  (1^  at  1.25) 
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Figure  9.  Double  Perturbation,  stable  (I  and  I  at  1.25). 
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Figure  11.  Triple  Perturbation,  stable  (I^,  I  and  M  at  1.25) . 


C.  STABILITY  WITH  d  AT  a. IB 

yx 

Similar  to  Table  IV,  Table  V  lists  the  initial  parcuneters 

for  investigation  of  the  model  with  d  at  0 . IQ  and  d 

yx  xy 

at  l.QO.  Graphically,  for  a  single  perturbation,  this  form 
of  the  model  produces  similar  results  as  before.  The  excep¬ 
tion  is  given  in  Figure  13  where  only  1^  is  perturbed 

(trial  3) .  M  and  M  can  be  seen  to  return  to  equilibrium 
X  y 

without  oscillation  as  had  been  indicated  in  the  previous 
investigation.  Again  all  trials  exhibit  the  dependent 
variables  damping  out  in  the  time  allowed  and  since  no 
significantly  differently  graphical  forms  were  produced  they 
are  not  included  in  this  section. 

D.  INSTABILITY  WITH  d  AT  1.4Q 

xy 

Table  VI  lists  the  initial  parameters  for  the  unstable 

case  where  d  =  1,4Q  and  d  is  held  at  l.QO.  The 
xy  yx 

additional  column  on  this  table  represents  which  force  term 

M  or  M  is  increasing  without  bound  while  the  other  term 
X  y 

is  decreasing.  It  should  be  noted  that  some  variables  were 
given  initial  values  of  l.QQl  to  get  comparative  plots  over 
time,  i.e.  using  1.25  as  an  initial  parauneter  in  these  cases 
caused  either  or  to  increase  or  decrease  severely 

immediately. 

Graphically  three  general  forms  of  instability  were 
demonstrated  where  only  one  dependent  variable  was  perturbed. 
These  are  presented  in  Figures  14,  15,  and  16  (trials  1 
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TABLE  V 


INITIAL  CONDITIONS  WITH 

d  AT  l.Q  AND  d  AT  Q.l  (STABLE) 
xy  yx 


Time  for 


Trial 

# 

Time 

Increment 

III 

90%  Damp 
Out  (Sec) 

1 

Q 

to 

50 

0.05 

1,25 

1.00 

1.00 

l.QO 

2.0 

2 

Q 

to 

50 

0.05 

1.00 

1.25 

1.00 

1.00 

21.0 

3 

Q 

to 

50 

0.05 

1.00 

1.00 

1.25 

1.00 

21.0 

4 

Q 

to 

50 

0.05 

1.00 

l.QO 

1.00 

1.25 

19.0 

5 

Q 

to 

50 

0.05 

1.25 

1.25 

1.00 

l.QO 

21.0 

6 

0 

to 

50 

0.05 

1.25 

1.00 

1.25 

1.00 

20.0 

7 

0 

to 

50 

0.05 

1.25 

1.00 

1.00 

1.25 

18.0 

8 

0 

to 

50 

0.05 

1.00 

1.25 

1.00 

1.25 

4.0 

9 

Q 

to 

50 

0.05 

1.00 

1.00 

1.25 

1.25 

34.0 

10 

0 

to 

50 

0.05 

1.00 

1.25 

1.25 

l.QO 

6.0 

11 

0 

to 

50 

0.05 

1.25 

1.25 

1.25 

l.OQ 

8.0 

12 

0 

to 

50 

0.05 

1.25 

1.25 

1.00 

1.25 

6.5 

13 

0 

to 

50 

0.05 

1.25 

1.00 

1.25 

1.25 

33.0 

14 

Q 

to 

50 

0.05 

l.OQ 

1.25 

1.25 

1.25 

12.5 
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TABLE  VI 


-i 


1 


INITIAL  CONDITIONS  WITH 

d  AT  1.4  AND  d  AT  1.0  CUNSTABLE) 
xy  yx 


Trial 

#  Time  Increment  '^'x 


1 

0 

to 

50 

in 

o 

• 

1.25 

2 

Q 

to 

50 

.05 

l.OQ 

3 

0 

to 

50 

.05 

l.QQ 

4 

0 

to 

50 

.05 

1.00 

5 

0 

to 

50 

.05 

1.25 

6 

0 

to 

50 

.05 

1.001 

7 

0 

to 

50 

.05 

l.QQl 

8 

0 

to 

50 

.05 

1.00 

9 

0 

to 

50 

.05 

l.OQ 

10 

0 

to 

50 

.05 

1.00 

11 

0 

to 

50 

.05 

1.25 

12 

0 

to 

50 

.05 

1.25 

13 

0 

to 

50 

.05 

1.001 

14 

0 

to 

50 

.05 

1.00 

Increasing 
Force  Term 

1.00 

1.00 

l.OQ 

1.25 

l.QQ 

1.00 

«x 

l.OQ 

l.QQl 

1.00 

“y 

l.QQ 

1.00 

1.001 

“y 

1.25 

1.00 

l.OQ 

"x 

l.OQ 

1.001 

l.QO 

“y 

l.QQ 

1.00 

1.001 

"y 

1.00 

1.001 

1.001 

«y 

l.QOl 

1.00 

1.001 

«x 

1.01 

1.01 

1.00 

“y 

1.25 

1.25 

1.00 

«x 

1.25 

1.00 

1.25 

«x 

1.000 

1.001 

1.001 

«y 

l.QOl 

1.001 

1.001 

I 
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Uii.  MU ippviipffllPIIHIPliV 


through.  41 .  The  correapondijig  phase  plots  are  given  in 
Figurea  17,  18,  19,  and  20.  Specific  values  for  double 
and  triple  perturbation  did  not  demonstrate  any  unusual 
graphic  forms  of  instability  and  generally  followed  those 
graphic  representations  displayed  for  single  perturbations. 

E.  INSTABILITY  WITH  d_  AT  1.4Q 

yx 

Table  VII  gives  the  initial  values  for  trial  runs  holding 

d  at  l.QO  and  d  at  1.4Q.  Similar  graphic  forms  were 
xy  yx 

produced  in  this  investigation  as  those  produced  in  the 
preceding  investigation  of  instability.  Two  exceptions  are 
given  in  Figures  21  and  22  (trial  7)  and  Figures  23  and  24 
(trial  10)  .  It  also  should  be  noted  that  in  this  case 
perturbing  and  1^  together  causes  to  increase. 

In  the  previous  case  perturbing  these  variables  together 
drove  upward,  an  interesting  result. 

F.  GENERAL  OBSERVATIONS 

The  cases  investigated  in  sections  B  and  C  of  this 

chapter  clearly  indicate  that  the  Dynamic  Model  of  Modern 

Military  Conflict  does  exhibit  stability  as  predicted  in 

Chapter  3.  Furthermore,  it  also  may  behave  in  an  unstable 

manner  as  indicated  in  sections  O  and  E  above.  A  closer 

look  at  the  unstable  cases  also  reveals  bifurcation  relative 

to  the  values  for  d  and  d  ,  initial  values  for  the 

xy  yx 

dependent  variables  and  the  resulting  increasing  term. 
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Figure  17.  Single  Perturbation,  unstable  phase  plot  (1^^  at  1.25). 


Figure  18.  Single  Perturbation,  unstable  phase  plot  (M  at  1.25) . 


7.96E-02  .310  .5«0  .770  1.00 


Figure  20.  Single  Perturbation,  unstable  phase  plot  (M  at  l.OQl). 


TABLE  VII 


INITIAL  CONDITIONS  WITH 


d 

xy 

AT  1.0 

and 

AT  1.4 

(UNSTABLE) 

Trial 

# 

Time 

Increment 

M 

X 

I 

M 

Increasing 
Force  Term 

1 

0 

to 

50 

.05 

1.25 

l.OQ 

1.00 

l.OQ 

M 

X 

2 

0 

to 

50 

.05 

1.00 

1.25 

1.00 

l.OQ 

M 

X 

3 

0 

to 

50 

.05 

1.00 

l.OQ 

1.25 

l.OQ 

M 

y 

4 

0 

to 

50 

.05 

l.OQ 

1.00 

1.00 

1.25 

M 

y 

5 

0 

to 

50 

.05 

1.25 

1.25 

1.00 

l.OQ 

M 

X 

6 

0 

to 

50 

.05 

1.01 

l.OQ 

1.01 

1.00 

M 

y 

7 

0 

to 

50 

.05 

1.01 

l.OQ 

l.OQ 

1.01 

M 

y 

8 

Q 

to 

50 

.05 

1.00 

l.OQ 

1.01 

1.01 

M 

y 

9 

0 

to 

50 

.05 

1.00 

1.01 

1.00 

1.01 

M 

X 

10 

0 

to 

50 

.05 

1.00 

1.01 

1.01 

1.00 

^x 

11 

0 

to 

50 

.05 

1.25 

1.25 

1.25 

1.00 

M 

X 

12 

0 

to 

50 

.05 

1.25 

1.25 

1.00 

1.25 

«x 

13 

0 

to 

50 

.05 

1.001 

l.OQ 

1.001 

1.001 

‘"y 

14 

0 

to 

50 

.05 

l.OQ 

1.001 

1.001 

1.001 

M 

V 
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Figure  21.  Double  Perturbation,  unstable  (I  and  M  at  1.01) 


Figure  22.  Double  Perturbation,  unstable  phase  plot  (I 


****»«  M  ••••*>  II  iixii«  ai 


Figure  23.  Double  Perturbation,  unstable  (M  and  I  at  1.01) 


lO't  I9’l  tO't 


Figure  24.  Double  Perturbation,  unstable  phase  plot  (M  and  I  at  1.01). 


A  simple  display  of  th±s  bifurcation  is  seen  when  only 
one  element  (dependent  variable),  is  perturbed  as  in  trials 
1  through  4  in  Table  VI  and  VII .  Relative  to  Table  VI ,  d 

xy 

is  set  at  1.4Q  auid  d  at  l.OQ.  This  gives  Y  an 

yx 

advantage  over  X  in  c3  effectiveness,  however,  it  can  be 

seen  that  in  trials  1  and  2  of  this  table,  X  can  overwhelm 

Y  (M^  increasing,  decreasing)  if  given  an  advantage  in 

either  information  (I  )  or  forces  (.M  )  .  Similarly,  relative 

to  Table  VII  with  d  set  at  l.OQ  and  d  at  1.4Q,  Y 

xy  yx 

can  be  seen  to  overwhelm  X  in  trials  3  and  4 .  Evidences  of 
bifurcation  are  a  bit  more  complicated  in  the  remaining  trials 
of  these  tables. 

In  Table  VI,  trials  5,  9,  11,  and  12  exhibit  this 

3 

bifurcation  where  y  has  an  advantage  in  C  but  is  over¬ 
whelmed  by  X.  Similarly  in  Table  VII,  trials  6,  7,  8,  13 
and  14  show  X  being  overwhelmed  by  Y,  even  though  it  has 
a  C  advantage.  Trail  10  is  not  a  bifurcation  example,  but 
as  mentioned  in  section  E,  does  produce  an  interesting  result. 
Specifically,  from  Table  VI,  a  slight  advantage  in 

3 

■and  ly  results  in  Y  overwhelming  X  when  Y  has  a  C 

advantage.  Whereas  from  Taible  VII,  the  slight  advauitage  in 

M  and  I  results  in  X  overwhelming  Y  when  X  has  the 
X  y 

c3  advantage.  This  seems  consistent  with  common  sense. 


V.  EQUILIBRIUM  ANALYSIS 


All  the  preceding  analysis  of  stability  was  predicated 
upon  a  constant  equilibrium  vector  of  1  for  ease  of  calcula¬ 
tion.  Another  question  that  should  be  answered  at  this  point 
is  "what  happens  to  equilibrium  values  if,  for  instance,  when 

replenish  rates  (Q  ,R  ,Q  ,R  )  are  kept  constant  and  d  and/or 

X  X  y  y  xy 

dy^  is  allowed  to  vary?"  To  answer  this  question  the  equi¬ 
librium  equations  and  initial  analysis  conditions  are  brought 
back  into  the  analysis. 


A.  METHOD  OF  ANALYSIS 

First,  the  original  equilibrium  equations  are  set  to  0  and 
solved  for  the  replenishment  rates.  This  result  is  then  used 


to  solve  the  same  set  of  equations  for  equilibrium  as  a  func¬ 
tion  of  d  and  d  . 

r  * 

B.  EQUILIBRIUM  RELATION  TO  d_  AND  d  ^ 

xy  yx 

Each  equilibrium  equation  is  set  to  0  and  solved  for 

replenishment  rates. 

CM  -al  -yM  -A+Q=0 
XX  xe  xy  ye  'xy  ye  x  '^x 

0(1)  -  (1)  (1)  -  0(1)  -  0.5  +  Q^  =0 

=  1.5  (Eqn  5.1) 


-dl  -(6+d)I  -(S  +b)M  -B+R»0 

XX  xe  xy  xy'  ye  '^xy  xy'  ye  x  x 

0(1)  -  (0  +  1  -(1  +  0)  1  -  0.5  +  R^  »  0 


"x  -  '^xy  " 


(Eqn  5.2) 


-a  I  -yM  +C  M  -A+Q=Q 
yx  xe  'yx  xe  yy  ye  y  y 

Q  (11  -  1  (11  +  Q  (.11  -  0.5  +Qy  «  0 

Qy  =  1-5 


(Eqn  5.3) 


-(6  +d)I  -(6  +b)M  -dl  -B+R«0 

'  yx  yx^  xe  '^yx  yx'  xe  yy  ye  y  y 

-  (0  +  dy^ll  -  (  1  +  0  11  -  0  (1)  -  0.5  +  Ry  »  0 


Ry  -  V  1-5 


(Eqn  5.4) 

Utilizing  these  results  the  same  equations  are  solved 


again  for  the  equilibrium  points  starting  with  some  arbitrary 

initial  values  for  d  and  d  /  d'  and  d'  respectively 

xy  yx  xy  yx 

to  fix  the  replenishment  rates. 

CM  -a  I  -yM  -A+Q=0 
XX  xe  xy  ye  'xy  ye  x  '^x 

0  (1)  (lyg)  -  0.5  *  0 

I,,  .  Q,  -  0.5 

Using  the  results  of  equation  5.1 

(Eqn  5.5) 


xy 


^ye 

=  1 

•d  I 

(5 

+  d  )I  ^  -  (3 

XX  xe 

'  xy 

xy  ye 

■(01I„  - 

(0  + 

J  ^  •  &«»»» 

'  xe 

xy  ye  xy 

M 

*  -d  I  -  0 . 

ye 

xy  ye 

r  the  result  of  equations  5.2 

■  1.5  for 

«x 

*  -d  -  Q.5  + 

ye 

xy 

M 

»  1.0  -  (d„  - 

ye 

xy 

xy  xy  ye 
+  0)  -  ( 


=  0 


'  xy 


(Eqn  5.6) 
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0 


I 

Y 

M  + 

C  M 

yx  xe 

'yx 

xe 

yy  ye 

"»Ixe 

-  (X)V  * 

e 

O 

=  Q 

-  Q.5 

xe 

y 

=  0 


Using  the  results  o£  equation  5.3 


“xe  =  ^ 


(Eqn  5.7) 


-‘«yx  *  -  '8yx  *  ‘=yx>"xe  '  '^yy^ye  '  \  ^  "y  =  “ 

-«>  *  V^xe  -  <  1  *  0  ).«xa  -  "»V  ■  “•^^’'y  '  “ 
Using  the  results  of  equation  5.7 


1.5  -  R 


xe 


yx 


Using  the  results  of  equation  5.4  and  substituting 
‘^yx*  ^ 


Adding  d  -  d  to  the  numerator 
yx  yx 

(d  -  d*  ) 

I  =  1  -  y;  yx' 

*®  V 


C.  ASSIMILATION 

It  suffices  to  say  then  that  I^^  and  M^^  are  not 

effected  by  changes  in  d  or  d  •  On  the  other  hand 

xy  yx 

M  varies  in  the  opposite  direction  to  deviations  from 
ye 

its  nominal  value.  I  varies  in  the  opposite  direction  to 

X6 

deviations  from  its  nominal  value  normalized  by  d  . 

xy 
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More  specifically,  equation  5.6  indicates  that  force  Y 
would  require  fewer  general  purpose  resources  to  hold  force 
X  to  1  provided,  force  Y  increased  its  enhamcement 
beyond  its  initial  value.  In  a  similar  fashion  equation  5.8 
indicates  force  X  would  not  require  as  much  information  to 
hold  force  Y  to  1  provided  force  X  increased  its 
enhancement  beyond  its  intial  value.  It  should  be  noted 
that  these  equations  may  not  be  true  in  general,  but  rather 
may  only  be  true  for  the  specific  initial  values  used  in 
this  analysis. 


1 


VI.  SUMMARY  AND  CONCLUSIONS 

This  thesis  has  described  the  initial  investigation  of 
the  Dynamic  Model  of  Modern  Military  Conflict.  Following  a 
brief  description  of  the  model,  investigation  of  stability 
was  undertaken  using  a  method  of  analysis  previously  prescribed 
for  the  study  of  deterministic  systems  such  as  ecosystems. 
Through  this  investigation  it  has  been  demonstrated  that  the 
model  does  in  fact  behave  as  one  might  expect,  given  a  set  of 
somewhat  reduced  initial  parameters.  In  this  respect,  this 
research  endeavor  did  not  attempt  to  fully  prove  or  disprove 
the  model  but  rather  its  purpose  was  to  permit  observation  of 
the  model  behavior  as  am  initial  step  toward  further  study. 

A.  GENERAL  OBSERVATIONS 

It  can  be  seen  that  the  Dynamic  Model  of  Modern  Military 

Conflict,  given  a  set  of  asymmetrical  initial  parameters 

exhibits  behavior  that  corresponds  to  a  status  quo  (stability) 

or  a  dominance  of  one  force,  combat  force  if  you  will,  over 

another  (instability), .  The  significance  of  this  observation 

lies  in  the  plausible  use  of  this  model  in  predicting  outcomes 

of  military  conflict.  Particularly,  the  incorporation  of 

3  3 

terms  relating  to  counter  C  and  C  enhancements  along  with 
more  tradition.  1  terms  such  as  loss  rates,  replenishment 
rates  and  Lanchester  coefficients,  adds  a  contemporary  flavor 
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to  this  model.  The  actual  employment  of  this  model  in 
playing  "what  if?"  type  wargames  however  cannot  be  undertcUcen 
as  yet  without  additional  research  including  the  remaining 
teinns  of  this  model. 

B.  IMPLICATIONS  FOR  FUTURE  RESEARCH. 

Specifically,  the  following  areas  for  additional  research 
are  suggested  based  on  the  initial  investigation  described  in 
this  thesis: 

1.  Location  of  other  equilibrium  points  besides  the 
vector  1  utilized  herein. 

2.  Determination  of  an  appropriate  Lyapunov  function 
suitable  for  predicting  global  stability. 

3.  Prediction  of  bifurcation  through  an  analytical 
approach. 

4.  Further  investigation  of  model  behavior  relative  to 
variations  of  other  initial  parameters  (i.e.  counter  C3 
coefficients )  . 
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